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Intermodulation atomic force microscopy (ImAFM) is a mode of dynamic atomic force microscopy 
that probes the nonlinear tip-surface force by measurement of the mixing of multiple tones in a 
frequency comb. A high Q cantilever resonance and a suitable drive comb will result in tip motion 
described by a narrow-band frequency comb. We show by a separation of time scales, that such 
motion is equivalent to rapid oscillations at the cantilever resonance with a slow amplitude and phase 
or frequency modulation. With this time domain perspective we analyze single oscillation cycles in 
ImAFM to extract the Fourier components of the tip-surface force that are in-phase with tip motion 
(Fj) and quadrature to the motion [Fq). Traditionally, these force components have been considered 
as a function of the static probe height only. Here we show that Fj and Fq actually depend on both 
static probe height and oscillation amplitude. We demonstrate on simulated data how to reconstruct 
the amplitude dependence of Fi and Fq from a single ImAFM measurement. Furthermore, we 
introduce ImAFM approach measurements with which we reconstruct the full amplitude and probe 
height dependence of the force components Fi and Fq , providing deeper insight into the tip-surface 
interaction. We demonstrate the capabilities of ImAFM approach measurements on a polystyrene 
polymer surface. 
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INTRODUCTION 

Since its invention T atomic force microscopy (AFM) has developed into one of the most versatile techniques in 
surface science. At length scales ranging from micrometers down to the level of single atoms, AFM-based techniques 
are used to image[^-0], measure[51[B] and manipulate mater[3{5] at an interface. As an imaging tool, the goal of AFM 
development has been increasing spatial resolution and minimizing the back action force from the probe on the sample 
surface. A major advancement in this regard was the development of dynamic AFM^IO in which a sharp tip at the 
free end of the AFM cantilever oscillates close to the sample surface as depicted in figure [l] In order to achieve stable 
oscillatory motion, an external drive force is applied to the cantilever which is usually purely sinusoidal in time with 
a frequency that is close to the resonance frequency of the first flexural eigenmode of the cantilever. The high quality 
factor of the resonance ensures that the responding motion of the tip is approximately sinusoidal in time, with the 
same frequency as the drive signal [TTl I12| . Such periodic motion is best analyzed in the frequency or Fourier domain, 
where the motion is well described by one complex-valued Fourier coefficient at the drive frequency. This motion 
has a corresponding Fourier coefficient of the tip-surface force, which can be expressed in terms of two real-valued 
components, Fj which is in-phase with the motion and Fq which is quadrature to the motion. At a fixed probe height 
h above the surface, the two force quadratures Fj and Fq give only qualitative insight into the interaction between 
the tip and the surface [T5] and most quantitative force reconstruction methods are based on a measurement of Fj 
and Fq at different h [TiHT^ . 

In order to increase the accessible information while imaging with AFM, a variety methods have been put forward 
where amplitude and phase at more than one frequency is analyzed. These multifrequency methods can be divided in 
to two general groups: Those using only Fourier components with frequencies close to a cantilever resonance, and those 
which use off-resonance components. Off-resonance techniques typically measure higher harmonics of the tip motion 
which allows for a reconstruction of time-dependent surface forces acting on the tip. Due to the lack of transfer gain 
off resonance, these off-resonance components have small signal-to-noise ratio and their measurement requires special 
cantilevers pOj. high interaction forces[21 or highly damped environments 22J. To increase the number of Fourier 
components with good signal-to-noise ratio, on-resonance techniques utilize multiple eigenmodes of the cantilever |231- 
126] . However, accurate calibration of higher cantilever modes remains complicated since additional knowledge about 
the cantilever is required. Both on and off-resonance techniques require broad-band detection of the cantilever motion, 
which implies a sacrifice in the sensitivity and gain of the motion detection system. 

To mitigate these problems, we have developed narrow-band intermodulation AFM (ImAFM) which analyzes the 
response only near the first flexural eigenmode. In general ImAFM utilizes frequency mixing due to the nonlinear tip- 
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surface interaction. A drive signal that comprises multiple frequency components is used for exciting the cantilever, 
which will exhibit response not only at the drive frequencies, but also at frequencies that are linear integer combinations 
of the drive frequencies, 

'^imp = miWi + ■m2U}2 + . . . + niMUJM (1) 

where mi , m2 , . . . mjv/ G Z and wi , a;2 , . . . , ujm are the drive frequencies. These new frequency components are called 
intermodulation products (IMPs) and one usually defines an order for each IMP which is given by \mi \ + \m2 \ + . . . + 
|mAf |. If all frequencies in a signal are integer multiples of a base frequency Auj the signal is called a frequency comb. 
The nonlinear tip-surface interaction maps a drive frequency comb to a response frequency comb, both having the 
same base frequency Aw. Different drive frequency combs can be used to place many response frequency components 
close to a resonance of the cantilever where they can be detected with good signal-to-noise ratio. In general the 
drive and response frequency combs could encompass more than one eigenmode of the cantilever. For a drive signal 
consisting of only two frequencies symmetrically placed around the first flexural resonance frequency as illustrated in 
figure [l] the response is concentrated in the narrow band around the first resonance for which well accurate calibration 
methods exist [27H^. 

In what follows, we will focus on this particular case which we call narrow-band ImAFM. However, we want to 
emphasize that drive schemes which generate response in more than one frequency band are also possible. We have 
previously shown how the individual amplitudes [^Hj and phases [301 of the IMPs in the narrow band around the first 
flexural resonance can be used for imaging. Furthermore, a polynomial reconstruction of the tip-surface forcef3ll 132] 
and a numerical fit of the parameters of a force mo del [55] are possible by analysis of the data in the frequency domain. 
Here, we consider the meaning of the narrow-band intermodulation response comb in the time domain, which leads 
to a physical interpretation of the intermodulation spectrum in terms of the in-phase force component Fi and the 
quadrature force component Fq. 



RESULTS AND DISCUSSION 
Time domain interpretation of narrow-band frequency comb 

Figure [2][a) portrays the amplitudes of the components of a narrow band frequency comb. While we are plotting 
only the amplitude at each component, it is understood that each component also has a phase. The frequency comb 
is characterized by the center frequency ui of the comb, a base frequency Alu and a finite number N of Fourier 
components at discrete frequencies. Without loss of generality, we assume N to be even. The center frequency ui can 
be described in terms of the ratio n = "/aw and the bandwidth of the comb is given by AA^ = A — 1. The A^ discrete 
frequencies in the band are represented by an integer frequency index n such that 

u)n = nAuj (2) 

where n takes consecutive integer values between n — ^^/2 and n + ^^/2. In the time domain the corresponding 
real- valued signal x{t) is then given by the Fourier series 

x{t)^ ^ i„e"^"*+ Ce"""^"* (3) 

n=ft-A«/2 n=n-AN/2 

where x„ are the complex Fourier components in the narrow frequency band and the star denotes complex conjugation. 
The center frequency w is usually much bigger than the base frequency Aw, 

to > Auj. (4) 

Therefore, the time domain signal x{t) exhibits two different time scales: a fast time scale Tfast = ^V'^ and a slow 
time scale Tgiow = ^tt/aw. To separate these two time scales we factor out a rapidly oscillating term at the frequency 
w from the Fourier series in equation ([s]), 

x{t) = Xne'^"'^"~")*e'"* Y Ke~"^"'^"~"'''e-''^\ (5) 

n=n— AJV/2 n=fi— AN/2 

ft + AJV/2 n+&N/2 

n=n-AN/2 n=fl-AN/2 
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FIG. 1: Sketch of the basic experimental setup in narrow-band ImAFM. In the absence of a drive signal the tip is at rest at 
the static probe height h. The spectrum of the drive comb consist of two frequency components spaced by Ao; and centered 
at the first fiexural resonance frequency cjq of the cantilever which is much higher than the comb base frequency Aw (here 
cjo = 600 Atj). The driven tip oscillates with amplitude A and interacts with the sample surface. The instantaneous tip position 
z is measured in the rest frame of the sample surface. The corresponding deflection signal is detected by an optical lever system 
and is concentrated to a narrow band around cjo, as the drive signal. In this band new frequency components spaced by Aoj 
are present which are generated by the nonlinear tip-sample interaction Fts- Outside the narrow band at u)o there is only little 
response in bands at integer multiples of ujq. 



Since N is even, 2n is an odd integer number and we can define a new sum index 

n' = 2n-2n 

which increase in steps of An' — 2 and the summation limits become 

n'l = 2 (n - ^) - 2n = -{N - 1), 



n'j^ = 2(n 



2n =N -I. 



(7) 



(8) 
(9) 



Since N is even and n' increases only in steps of two, n' can only take odd values. Additionally, we define new Fourier 
coefficients 



such that the signal Fourier series becomes 



x{t) 



^n' — ^i(ri'+2fi) 



(10) 
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FIG. 2: The amplitude spectrum of a narrow band signal as a function of the Fourier index(a) is characterized by a finite 
number of Fourier components in a frequency band around the center frequency ijj — hAuj in a bandwidth that is given by the 
integer number AN. The spectrum of the corresponding time-dependent envelope function(b) is obtained by down-shifting the 
original spectrum in frequency space such that the shifted center frequency is zero. 



We identify the terms in parentheses as the Fourier series of a complex- valued time-dependent envelope function E{t) 
expanded in the base frequency ^'^/2, 



N-l 



and write the original signal x{t) as 



n' = -{N - 1) 
n' odd 



x{t) = E{t)e''^' +E*{t)e-' 



(12) 



(13) 



The envelope function E{t) was obtained by down-shifting the narrow intermodulation frequency band to a center 
frequency of zero (see figure [2|. If the maximum frequency in the Fourier seri es o f E{t) is much smaller than w, 
the envelope function E{t) varies slowly compared to the term e'"* in equation (13). When we represent E{t) by a 
time-dependent amplitude A(t) and a time-dependent phase (f)(t) such that 

E{t) ^ A{t)e"'''^*\ (14) 

the signal x(t) is completely described by a modulated oscillation amplitude and a modulated oscillation phase, 

x{t) = A{t)e"^'-'^e"^* + A{t)e-'^^'^e-'^' ^2A{t) cos {Cot + ,p{t)). (15) 

We would like to emphasize that the narrow-band frequency comb can also describe amplitude and frequency 
modulated signals. For frequency modulation we define an instantaneous oscillation phase 



e{t) ^ujt + (t){t) 



and an instantaneous oscillation frequency 



, , de 

u.{t)^-=u 



dt 



The instantaneous frequency shift 5uj compared to uj is then simply 



(16) 
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FIG. 3: The amplitude spectrum of an amplitude and frequency modulated signal(a). The response is concentrated in a narrow- 
band frequency comb with a center frequency much higher than the comb base frequency. In the time domain(b) the signal 
rapidly oscillates on a short time scale and the slow amplitude modulation is clearly visible. The time-dependent amplitude(c) 
and frequency(d) reconstructed from the envelope function are in excellent agreement with the actual modulation used for the 
signal generation in the time domain. 



In a small region around the time t' the signal x{t) can be obtained by a Taylor expansion to first order of the 
instantaneous phase 9 

x{t) = A{t) cos {e{t)) « A{t) cos ((w + 5u:{t')) t + {4>{t') - Su{t')t')) . (19) 

Thus, a narrow-band frequency comb can describe signals with frequency shifts that are periodic in time. 

To illustrate the complete description of a narrow-band signal by its envelope function, figure [3] shows the spectrum 
of an artificially constructed signal with sinusoidally modulated amplitude and frequency. Typical parameters from 
AFM experiments have been chosen for the amplitude and frequency modulation. The spectrum of the signal shown 
in figure [3] shows significant amplitudes at discrete frequencies in only a narrow band around 300 kHz. We down-shift 
the spectrum to determine the slowly varying envelope function from which we compute the time-dependent oscillation 
amplitude and frequency, both of which are in excellent agreement with the actual amplitudes and frequencies used 
for the signal generation. 
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FIG. 4: Sketch of a narrow band signal. On the slow time scale Tsiow the tip motion shows an amplitude modulation. On the 
fast time scale Tfast the signal rapidly oscillates. During each oscillation cycle the tip interacts with the surface for the time 
Tintcr during which the oscillation amplitude phase are approximately constant. 



To summarize, we have introduced a time domain interpretation of narrow-band frequency combs. If the center 
frequency of the band is much higher than the base frequency of the comb, we can separate a fast and a slow time 
scale in the time domain. On the fast time scale the signal rapidly oscillates at the center frequency. The slow time 
scale dynamics is given by the down-shifted intermodulation spectrum which describes a slow amplitude modulation 
and a slow phase or frequency modulation of the signal in the time domain. 



Physical interpretation of the tip motion and force envelope functions in ImAFM 



In ImAFM the measured frequency comb corresponds to a vertical motion z{t) of the tip which undergoes rapid 
oscillations at frequency ui with slowly varying amplitude and phase, 



z{t) = A{t) cos {uit + (l){t)) + h 



(20) 



where h is the static probe height above the surface and the amplitude A{t) and the phase (j){t) are determined from 
the complex- valued motion envelope function Ez{t) as 



m 



arg (E,{t) 



(21) 
(22) 



The envelope function E^it) was obtained directly from the measured motion spectrum using equation (12). 

Knowledge of the cantilever transfer function G and the applied drive force allows for converting the measured 
motion spectrum into the spectrum of the time-dependent tip-surface force acting on the tip. However, the force 
spectrum is incomplete since higher frequency components of the force are filtered out from the motion spectrum 
by the sharply peaked cantilever transfer function. The time-dependence of the corresponding partial force signal is 
described by the force envelope function Ep{t)^ 



F, 



partial 



(t) = EF{t)e''^' + E*p{t)e 



(23) 



where Ep(t) is determined by applying equation (12 1 to the partial force spectrum. To understand the physical 



meaning of the partial force we analyze the signals at the level of single rapid oscillation cycles in the time domain. 
During each oscillation cycle the tip interacts with the sample surface. This interaction is very localized (a few 
nanometers above the surface) compared to the oscillation amplitude (tens of nanometers) and thus the interaction 
time Tintor is short compared to the fast oscillation period Tfast which itself is much shorter than the period of the 
beat Tsiow (see figure |4|, 



Tintcr < Tfast ^ ?sl 



(24) 
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Therefore, amphtude A{t), phase 4>{t) and force envelope function Epit) can be considered to be constant during each 
interaction cycle and the motion and the partial force during the i-th tip oscillation cycle are given by 

= yl(^) cos (cj< + + (25) 

-^partial (^) - ^pC + E p € (25) 

where A'^*^ and e'P are constant and are determined at the time t'-*-' of the ith lower turning point of the tip 
motion 



. (^^*^) I (27) 
« = </)(tW) =arg(4(<W)) (28) 



Ef (t**') (29) 



The complete time-dependent tip-surface force during an interaction cycle is a force pulse which can be written as a 
Fourier series in the oscillation frequency Co as 

CXD 

j^i:Ute(*) = E (30) 



where the complex Fourier components Fn^ fulfill the relation f'^\ — Fn^* . Comparison of equation (|26|) with 



equation (30) yields 

(iW) (31) 



) = E^;^ = Ef 



which reveals that the first Fourier component F^^'^ of the force pulse during the ith oscillation cycle is given by the 
force envelope function Ep determined from the partial force spectrum. Since each lower motion turning point is 
associated with a unique amplitude ^'•^^ we can consider Fi as a function of the continuous variable A, 

f['^ ^ FiiA^^y) ^ Fi ^ Fi{A). (32) 

The amplitude-dependence of -Fi can then be uncovered by the analysis of all oscillation cycles during the time Tsiow 
To better relate motion and the force we compute the components of Fi{A) that are in phase with the motion (Fj) 
and quadrature to the motion (Fq). For a tip-surface force that only depends on the instantaneous tip position and 
velocity, 

Ft3 = Ft,(z,i), (33) 

we approximate the tip motion to be purely sinusoidal at frequency ui with amplitude A and without an additional 
phase. At fixed probe height h, the components Fj and Fq are given by two integral equations 

Fi{A,h) = - — / Fts{Acos{ijt) + h,-uiAsin{iDt))cos{ujt)dt (34) 



ast JO 



1 



FQ{A,h) = / Fts (Acos(wi) -I- /i, -wAsin((Dt)) sin(w<)dt (35) 

Tfast JO 

With these assumptions Fj becomes the so-called virial of the tip-surface force which is only affected by the conservative 
part of the tip-surface interaction[31] whereas Fq is related to the energy dissipated by the tip-surface interaction|llj. 
We note that through their dependence on tip position z and velocity z, the force components Fi and Fq are functions 
of both probe height h and oscillation amplitude A . However, usually they are considered as functions of the probe 
height h only. 

For an ImAFM measurement at fixed probe height, the amplitude-dependence of Fj and Fq can readily be obtained 
by defining a new force envelope function that is phase-shifted with respect to the motion by the angle 0, 

E'p{t) = FF(i)e""^^*) (36) 
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which we evaluate as real and imaginary parts at the times of the lower turning points of the motion 

FK^W) = Re(^^(iW)) , (37) 

i^Q(AW) = Im(i^(tW)) (38) 

With this interpretation of an intermodulation spectrum we are able to reconstruct the amplitude-dependence of the 
force quadratures Fj and Fq which are independent of details of the tip motion on the slow time scale. Due to this 
independence, Fj and Fq are the input quantities for nearly all force spectroscopy methods in dynamic AFM and 
thereby they form the basis of quantitative dynamic AFM. 



Force quadrature reconstruction from simulated data 

To demonstrate the accuracy of the Fj (A) and Fq {A) reconstruction from ImAFM data we simulate the tip motion 
in a model force field. We excite the tip with two frequencies close to the first flexural resonance frequency of 
cjQ = Stt • 300 kHz which allows us to model the cantilever as a single eigenmode system for which the tip dynamics 
are described by an effective harmonic oscillator eauation |351 136| 



z + + ujI{z - h) = {Fi cos(wit) + Fa cos(w2t) + Fts(z, i)) 



(39) 



where Q = 400 is the quality factor of the resonance, kc = 40 N/m is the mode stiffness and h — 20 nm is the 
static probe height above the surface. The drive strengths Fi and F2 at the frequencies wi = 2n ■ 299.75 fcHz 
and a;2 = 27r • 300.25 kKz are chosen such that in the absence of a tip-surface force, the tip oscillation amplitude is 
sinusoidally modulated between nm and 30 nm . For the tip-surface force F^g we assume a van-der-Waals-Derjaguin- 
MuUer-Toporov (vdW-DMT) force with additional exponential damping as which is defined as 



Fts{z,z) 



HR 
'6(ao-|-2:)2 
fffl I 4 

Qai 3 



7 exp (—2/2^) i z > 

^£:*(-z)3/2 -7cxp(-2/2-,)^ z<Q 



(40) 



where H = 2.96 • 10~^ J is the Hamaker constant, i? = 10 nm is the tip radius, 7 — 2.2 • 10~^ Ns/m is the damping 
constant, = 1.5 nm is the damping decay length and E* = 2.0 GPa is the effective stiffness. For the numerical 
integration of equation (39 1 we use the adaptive step-size integrator cvode[37] with root detection to properly treat 
the piecewise definition of the tip-surface force in equation (40). From the simulated tip motion we deterrnine the 



motion and the force envelope functions Ez and Ep and reconstruct Fj{A) and Fq{A) according to equations (37) and 
(38). As shown in figure ([s]) the reconstructed curves are in excellent agreement with the curves directly computed 



with equations (34) and (35) from the model force used in the simulations. 



Probing the force quadratures 

The force quadratures Fi and Fq are the basic input quantities for a variety of force reconstruction techniques |14l- 
[TU] . Over the last decade the dominate paradigm was to consider F/ and Fq as functions of the static probe height h 
only, and only one oscillation amplitude A was considered at each probe height. F/ and Fq are however functions of 
both h and A as seen in the two-dimensional color maps shown in figure |6] for the vdw-DMT force with exponential 
damping used in the previous section. In order to emphasize the interaction region near the point of contact, data in 
the h-A plane with F/ < —8 nN are masked with white. 

In both frequency modulation AFM (FM-AFM) and amplitude modulation AFM (AM- AFM) F/ and Fq are usually 
probed by a slow variation of the probe height h with fixed oscillation amplitude at each height. To measure Fj and 
Fq in FM-AFM the oscillation frequency shift and the drive force are recorded as the static probe height is slowly 
varied (frequency-shift-distance curves). Active feedback is used to adjust both the drive power and drive frequency, 
as to keep the response amplitude and phase constant. The obtained frequency shifts and drive forces can then be 
converted into the force quadratures |38( I39j so that the measurement corresponds to a measurement of F/ and Fq 
along a path parallel to the /i-axis in the /i-A-plane (see figure [6]) . 

In AM- AFM the oscillation amplitude and phase with respect to the drive force are measured as a function of the 
static probe height h (amplitude-phase-distance curves) and are then converted into values of the force quadratures}!^] . 
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FIG. 5: The Fi{A) and Fq{A) curves reconstructed from simulated tip motion in ImAFM. The reconstructed curves are in 
good agreement with the actual curves directly determined from the model force used in the simulations. 

In contrast to FM-AFM, the oscillation amplitude is free to change during the measurement and thus the AM-AFM 
measurement path in the h-A-plane is more complicated. The path shown in figure |6] was obtained by simulating 
the AM-AFM tip dynamics with cvode. In the simulations we used the same cantilever and force parameters as in 
the previous section with a drive signal at only one frequency of uJd = 300 kHz. As is often the case with AM-AFM, 
the amplitude-phase-distance curve exhibits an abrupt amplitude jump due to the existence of multiple oscillation 
states [10]. This instability is frequently observed in experiments, and it makes the reconstruction of tip-surface forces 
rather difficult. 

In contrast to FM-AFM and AM-AFM, ImAFM allows for a measurement of Fj and Fq at fixed static probe 
height, along a straight path parallel to the A-axis in the h-A plane as shown in figure [6] for the simulation of the 
previous section. Each of these three measurement techniques probes the tip-surface interaction along a different path 
in the h-A plane. With ImAFM however, the measurement can be rapidly preformed in each point of an image, while 
scanning with normal speed |26j allowing for unprecedented ability to analyse the tip-surface force while imaging. The 
ImAFM spectral data, which is concentrated to a narrow band near resonance, is a complete representation of the 
measurable tip motion, because there is only noise outside this narrow frequency band. Thus the method optimally 
extracts the signal for compact storage and further analysis. 

We note that the ImAFM path provides an equivalent amount of information as frequency-shift-distance or 
amplitude-phase-distance curves. This implies that for a single scan ImAFM image information equivalent to a 
frequency-shift-distance curve or an amplitude-phase-distance curve is available in every image point. Moreover, the 
ImAFM measurement does not suffer from amplitude jumps since the stiffness of the cantilever resonance prevents 
big amplitude changes from one single oscillation cycle to the next single oscillation in the beat tip motion. 

ImAFM approach measurements 

It is possible to acquire maps of Fj and Fq in the full h-A plane with a protocol we call ImAFM approach 
measurements. Similar to the measurement of frequency-shift or amplitude-phase-curves, the static probe height 
above the surface is varied by slowly extending the z-piezo toward the surface. However, in contrast to FM-AFM and 
AM-AFM measurements, the oscillation amplitude is rapidly modulated as the probe slowly approaches the surface. 
Because the height variation is much slower (order of seconds) than the amplitude modulation (order of milliseconds) , 
the probe height can be considered to be constant during each amplitude modulation. In this case each amplitude 
modulation reveals the amplitude dependence of Fj and Fq at a constant probe height. From the different probe 
heights Fj and Fq can be reconstructed can be reconstructed in the full h-A plane. With FM-AFM or AM-AFM 
such a measurement would require much longer measurement time since multiple surface approaches with different 
amplitudes would be required. With ImAFM all the data is acquired during a single surface approach. 

We use ImAFM approach curves to reconstruct Fj and Fq maps on a Polystyrene (PS) polymer surface. We 
perform a slow surface approach and from the acquired data we reconstruct the Fj and Fq maps shown in figure IT] 



10 




Probe height [nm] Probe height [nm] 



FIG. 6: Model Fi(h,A) and FQ{h,A) maps for the vdW-DMT force with exponential damping introduced in equation (40 1. 
The displayed measurement paths correspond to a frequency-shift-distance curve in FM-AFM, an amplitude-phase-distance 
curve in AM-AFM and an ImAFM measurement. In contrast to FM-AFM and AM-AFM, the static probe height is constant 
during an ImAFM measurement and the h-A plane is explored along a path parallel to the A axis. One should also note the 
amplitude jump along the AM-AFM path at a probe height oi h — 30 nm. 
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FIG. 7: Fj{h,A) and FQ{h,A) maps reconstructed from an ImAFM approach measurements on a PS surface. The z-piezo 
extension corresponds to a relative change of the probe height h above the sample surface. 



On the ft-axis we show the piezo extension since the absolute probe height cannot be defined unambiguously in an 
experiment. The areas in the h-A plane that were not explored are displayed as white areas. 

The boundary to the white area in the upper part of the plots represents the maximum oscillation achievable for 
the fixed drive power. The boundary to the white area in the lower part of the plot corresponds to the minimum 
oscillation amplitude during one modulation period (one beat). This lower boundary shows interesting variations 
with the piezo extension. Further away from the surface at larger piezo extension, only positive values of Fj are 
achieved, which corresponds to the tip oscillating in a region where the net conservative force is purely attractive. As 
the surface is further approached, Fj also takes on negative values, when the net force becomes repulsive. At probe 
heights between 7 and 13 nm the strongest repulsive force is experienced. Around h =13 nm and below, the attractive 
region vanishes, which may be the result of a change in the cantilever dynamics due to relatively long interaction time, 
or a change in the hydrodynamic damping forces due to the surrounding air close the sample surface. One should 
also note that at this piezo extension the minimum oscillation amplitude begins to increase again. A possible artifact 
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of the measurement method may result in this low amplitude region, if the motion spectrum is no longer confined to 
a narrow frequency band, as assumed in the analysis. 

The map of Fq characterizes the dissipative interaction between tip and sample. The dissipative tip-surface force 
can be much more complex than the conservative part of the interaction since dissipative forces do not only depend 
on the instantaneous tip position as with the conservative force. The Fq map can provide detailed insight into the 
nature of the dissipative interaction since the full dependence of Fq on probe height h and oscillation amplitude A is 
measured. The Fq map shows a lower level of force than the Fj map and it therefore appears more noisy, because the 
tip-surface interaction is predominately conservative. The small positive values of Fq which occur far from the surface 
would imply that the tip gained energy from the interaction with the surface. This may be an artifact, but another 
possible explanation is some sort of hydrodynamic mode above the surfaceHT]. In both the attractive and repulsive 
region of the in-phase force Fj , the quadrature force Fq is predominantly negative and it decreases as the surface is 
indented, corresponding to a increasingly dissipative tip-surface interaction. However, the maximum dissipation does 
not coincide with the maximum repulsive conservative force, and the energy dissipation is largest at peak amplitude 
for piezo extensions between 2 and 6 nm. Another interesting feature of the Fj map is the fine structure in the 
contact region, between h — 10 and 20 nm piezo extension. These small step-like changes of conservative force are 
not present in the smooth force model function, and could be indication that the dissipative forces are resulting in 
small, irreversible modifications of the sample surface. 

CONCLUSION 

We presented a physical interpretation of tip motion when described by a narrow-band frequency comb in ImAFM. 
We showed by separation of time scales that the time domain signal of a narrow-band frequency comb is completely 
characterized by a complex-valued envelope function and a rapidly oscillating term. The application of this time 
domain picture to ImAFM allows for the reconstruction of two force quadratures Fj and Fq as functions of the 
oscillation amplitude A. The quantities Fj and Fq can be considered as two-dimensional functions, depending on the 
both probe height and the oscillation amplitude. Within this framework we find a connection between frequency-shift- 
distance curves in FM-AFM, amplitude-phase-distance-curves in AM-AFM, and ImAFM measurements. Moreover, 
we introduced ImAFM approach measurements which allow for a rapid and complete reconstruction of Fj and -Fq 
in the full h-A plane, providing detailed insight into the interaction between tip and surface. We demonstrated the 
reconstruction of Fj and Fq maps experimentally on a PS polymer surface. We hope that the physical interpretation 
of narrow-band dynamic AFM presented here, will inspire new force spectroscopy methods in the future which take 
advantage of the high signal-to-noise ratio and the high acquisition speed of ImAFM. 

EXPERIMENTAL 

The PS sample was spin-cast from Toluene solution on a silicon oxide substrate. Both PS (M^, — 280 kDa) and 
Toluene were obtained from Sigma- Aldrich and used as purchased. The measurements where performed with a Veeco 
Multimode II and a Budget Sensor BS300A1-G cantilever with a resonance frequency of /o = 311.838 kHz, a quality 
factor of Q = 539.9 and a stiffness of = 29.5 N/m which was determined by thermal calibration [29^. We choose 
the two drive frequencies fi = 311.585 kHz and /2 = 312.085 kHz symmetrically around the resonance frequency and 
the drive strengths such that the free oscillation amplitude is modulated between 0.0 and 29.7 nm. The probe height 
is changed with a speed of 5.0 nm/s. 
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